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ABSTRACT 



In the recent literature on AGNs it has been often reported 
that spectra of Seyfert 1 Galaxies show resonant absorp- 
tion lines of Fe K which are redshifted from the rest frame 
position. Such lines are often found with marginal signif- 
icance but, if real, could potentially open up new avenues 
to study the circumnuclear gas in the black hole environ- 
ment. It is also extremely important to take them into ac- 
count in X-ray spectral analysis because of the influence 
they have in the correct estimation of spectral parame- 
ters, Fe Ka line in primis. An XMM-Newton observation 
of Mrk 335 is reported here as a case study: a narrow 
feature has been detected at 5.9 keV, i.e. with a redshift 
corresponding to a velocity of V'--^ 0.15c. Preliminary re- 
sults on the statistical significance of narrow absorption 
and emission lines in a sample of PG QSOs observed by 
XMM-Newton are also included. 
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1. INTRODUCTION 



Since Active Galactic Nuclei have been discovered, it has 
been postulated that the powering mechanism is likely to 
be the release of gravitational energy of matter accreted 
on a supermassive black hole. Observational evidence for 
this was found in the redshifted and broad Fe K disc line 
seen in bright Seyfert galaxies (1). The picture on hard 
X-ray spectra has been made even more complex by the 
recent detection of narrow lines shifted from their rest- 
frame position in the Fe K band spectra of many Seyfert 
1 galaxies which may affect modeling of the Fe K emis- 
sion line. Previous cases of highly ionised redshifted ab- 
sorption Fe K lines were in fact found superimposed to 
the broad wing of the Fe Ka line in NGC 3516 (2) with 
ASCA data, and in E 1821H-643 (H) with Chandra data. 
(I?) reported on the presence of an unshifted Fe K ab- 
sorption line superimposed on the relativistic Fe K line in 
IRAS13349H-2438, with XMM-Newton data. High confi- 
dence detections of such features would be of crucial im- 



portance in testing the black hole paradigm for AGN and 
would provide a new additional tool to be used alongside 
the broad Fe Ka line. In fact, although the exact nature 
of the energy shift of such lines is as yet unclear, the most 
likely scenario for producing the observed features would 
involve a combination of gas orbiting in highly relativis- 
tic motion and/or gravitational shifts of the photons ( 2; 5). 
With the advent of XMM-Newton and Chandra, the num- 
ber of absorption features in active galaxies spectra have 
considerably increased ('3, ^ Narrow energy- 
shifted emission lines have also been detected in the hard 
X-ray spectra of many AGN (9; 10; 11; l^|ll[l3). The- 
oretical models have predicted the possibility that Fe K 
emission lines from the disc can be observed with a nar- 
row profile if the X-ray reflection arises as a result of 
magnetic flares in locahzed regions on the disc Cl5.;.16.) . 



2. MRK 335 CASE STUDY 



2.1. Spectral analysis 

Mrk 335 is a bright Seyfert 1 Galaxy at z=0.026, which 
was observed by Xmm-Newton for about 30 ks. A pre- 
vious analysis was reported by ( 17), who found evidence 
for a broad Fe Ka line associated to an ionised reflection 
component. Here, the data from the pn camera are pre- 
sented. A fit on the 2-10 keV data with a simple power 
law model yields a steep spectrum, with r=2.13lQ S^. 
The pn residuals from 3-9 keV are plotted in Fig. [3 
the energy band of the Fe Ka emission line is pictured 
on the data showing the presence of broad excess in flux 
not only above the position of the neutral line (6.4 keV), 
but even up to ^ 7.3 keV and down to ~ 5.9 keV. A 
deficit of counts in a notch-shape is also present at ^ 
5.9 keV. A Gaussian line is added to the power law, with 
energy, width and flux free to vary. The line is highly 
significant with Ax^=33 for 3 d.o.f, indicating an ion- 
ized and broad line. The line parameters are found to 
be E=6.62,ta\l keV (rest-frame), o-=0.40t^;?3 keV and 
EW=245!ti26 Although the residuals shape may sug- 
gest the presence of another Gaussian line, any attempt to 
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Figure 1. Data to model ratio: the 2-10 keV spectrum 



is fitted by a power law with r= 2.15Z_q q^ . The plot is 
in the source rest frame and the Fe Ka energy band is 
labelled for clarity. 



fit the data with 2 emission lines failed. To fit the notch- 
shaped feature another Gaussian line with negative inten- 
sity has been added. The fit yields x^=738/731 d.o.f. and 
the lines parameters are found to be E=6.31^q'2o keV, 
a=Q.l%t°^ll keV, EW=468t?^g eV for the broad compo- 
nent and as for the narrow one, E=5.92to 05 keV with an 
EW=52 l]^g eV (measured in absorption with negative in- 
tensity with respect to the continuum). The width of the 
absorption line is unresolved with CCD resolution and 
therefore it is kept fixed to 50 eV. The improvement in 
is for 2 degrees of freedom, corresponding 

to a level of confidence higher than 99.7 percent accord- 
ing to the F-test. This is an extremely basic parametriza- 
tion of the spectrum, meant purely to show the main fea- 
tures in the spectral curvature above ^ 5 keV. The line 
profile in fig. appears complex not only for the pres- 
ence of the absorption feature, but also because the resid- 
uals show a double-peak structure. As reported before, 
fitting two emission lines is not required by the fit so 
we have included only one broad Gaussian in our basic 
parametrization. When the spectrum is fitted with two 
Gaussian lines (emission and absorption), the energy of 
the broad line is consistent with 6.4 keV. A close look to 
figE reveals that the profile is very different from a Gaus- 
sian and that in this case the use of such model could be 
quite misleading. The profile is asymmetrical and skewed 
suggesting that if there is a broad line it could be modi- 
fied by relativistic effects. We used the DISKLINE model 
in Schwarzschild metric ( 18) where the line parameters 
other than the energy are: q, the line emissivity index, 
where the line emissivity j is a function of the emission 
radius r according to j ocr"'?; the inner radius rin and the 
outer radius rout of the accretion disc which define the 
area of the disc where the line is emitted; the inclination 
of the disc /, set as to be the angle between the line of 
sight and the normal to the disc. Fitting the broad line in 
this way yields E=7.12t°J} keV and EW=407;i°^ 

^=3.98^o;7^, i=2l°t% and x^/d.o.f.=735/730. The ab- 
sorption line is included in this fit as a negative Gaussian, 



as previously described. 

The presence of a diskline suggests that a reflection com- 
ponent should be included. Since the line energy is 
clearly indicative of a high ionisation state for Iron, the 
XION model developed by (J9i) was used to fit the spec- 
trum. In this way, the reflected spectrum is computed in 
hydrostatic balance, taking into account the ionization in- 
stability in the disc. After choosing one of the available 
geometries, this model calculates the distance between 
the disc and the source of X-ray photons, the accretion 
rate, the luminosity of the X-ray source, the inner and 
outer disc radii and the spectral index as free parame- 
ters. Relativistic smearing is included for a non-spinning 
black hole. We resolved to assume the simplest geometri- 
cal configuration (lamppost). After adding an absorption 
Gaussian line, the model provides a fairly good fit, yield- 
ing x^/d.o.f. = 738/730. 

We try to model the absorption line by adding an appro- 
priate warm absorber to the best fitting reflection model 
XION. In order to do that, a grid of XSTAR photoion- 
ization model was generated with solar elemental abun- 
dances and turbulence velocity of 100 km/s. Then, such 
model has been incorporated in XSPEC as a table model 
with 3 additional free parameters: i) the column density 
Nu, ii) the ionization parameter ^, which describes the 
state of the photoionized medium; iii) a redshift parame- 
ter which includes all the redshift contributions, namely 
the cosmological redshift {Zsource), the velocity of the ab- 
sorber {zinfiow) and the gravitational shift which the gas 
may be subjected to, if close enough to the black hole 
iZgrav)- The best-fitting parameters are consistent with 
a very ionized absorber, with log^ ^ 3.9 ergs cm s^^ 
where it is most likely that only the H-like and He-like Fe 
ions survive. The fit yields x^/d.o.f.=735/728. Four main 
absorption features are imprinted on the continuum (see 
Fig. □: the Ka and K/3 transitions of Fe XXV and XXVI 
produce the absorption lines, but only the Ka ones are 
sufficiently strong to be interesting for our purposes here. 
The absorption line detected in the data is consistent with 
the Fe XXVI Ka, whereas the other ones are too weak to 
be detected at the CCD resolution. Most striking is that, if 
such feature is identified as we propose, it requires to be 
redshifted corresponding to an inflow velocity of ~ 0.14 
c. Such value is inferred by measuring the energy shift of 
the Gaussian line peak (~5.92 keV) from the rest-frame 
position of 6.90 keV. 



2.2. A simple model for the inflow 

To further investigate the inflow hypothesis, a simple 
physical model has been developed and used to synthe- 
sise X-ray spectra for qualitative comparison with the ob- 
served absorption feature (a more thorough description 
of the model is included in the paper Longinotti et al. in 
prep.) Spectra were synthesised for the 2-10 keV re- 
gion using a Monte Carlo radiative transfer code based 
closely on that discussed by (20). For the computation 
presented here, the power-law photon index was fixed at 
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Figure 2. Rest frame plot of the model used to fit the data 
on the left. The combination of the effects of the reflection 
spectrum and the highly photoionized gas is visible: the 
first one reproduces the spectral curvature between 5 and 
7 keV, the latter imprints redshifted absorption features 
on the continuum. In our data, we observe the Fe XXVI 
Ka only, the other ones being too weak. 

r = 2.2 (based on fits to the data) and the normahsa- 
tion fixed by requiring the total 2-10 keV X-ray flux to 
match the observed luminosity of 1.8 x lO"*"^ ergs s^^. 
The disc luminosity was set to 60 per cent of the Ed- 
dington luminosity for the central black-hole, close to the 
value L/Ledd = 0.62 reported for Mrk 335 by The 
soft-excess was modelled as a black-body with tempera- 
ture 1.3 X 10^ K and normalisation fixed to the power- 
law component. In contrast to that presented by ( :20l) . 
the code used here includes the full special relativistic 
expression for the Doppler shift and approximately ac- 
counts for gravitational redshift using 



^iv{r)/c)v = v'^l-2Rg/r (1) 

where v' is the frequency of a photon at radius r as mea- 
sured in the comoving frame, /i is the usual direction 
cosine, 7 = (1 — j and v is the photon fre- 
quency that would be recorded by an infinitely distant 
observer at rest relative to the black-hole. Other rela- 
tivistic effects are neglected. A simple spherical geom- 
etry for the gas in a radial inflow has been considered. It 
is assumed that the gas occupies a region which extends 
from inner radius ri,, to outer radius Tout from the central 
black-hole. Here only the main results from such model 
are described. From a qualitative comparison to the com- 
puted spectra, the data can rule out an inflow model of a 
continuous flow extending from the vicinity of the cen- 
tral black-hole to much greater distances. The left panel 
in Figure |3l shows the 3-9 keV spectrum computed for a 
model of this sort with ri„ = 20 Rg, rout = 10'^ Rg, and 
a mass infall rate $ = 0.2 Mq yr~^. This model predicts 
very few spectral features since the ionisation state of the 
gas is very high. The dominant feature is a broad inverse 
P Cygni Fe XXVI Ka line, but a weaker feature due to the 
K/3 line of the same ion appears at harder energies. The 
computed absorption line is very broad, extending from 



around 5.5 keV to ~ 6.5 keV, where the deepest absorp- 
tion occurs. The large line width in the model is therefore 
the result of the large radial extent of the absorbing gas 
and of the large velocity range present in the flow. If the 
infalling gas is instead distributed over a narrow range of 
radii, the model is found to be a good description of the 
data (right panel in Figure |3}- A possible configuration 
could be described as a discontinuous infalling of shells 
or blobs of gas at different densities. With the model 
considered here, this was obtained by limiting the radial 
range occupied by the flow. 



2.3. Significance of tlie narrow line in Mrk 335 

The problem of the reality and significance of narrow 
features such as the one detected in Mrk 335, has been 
pointed out in the most recent cases of narrow line de- 
tections by (3) for an absorption line, and by (13) for 
an emission line. These authors have employed realis- 
tic Montecarlo simulations for testing the reality of the 
lines, which have been found significant at a level in be- 
tween 2-3 cr. Moreover, the employment of the F-test 
to test the significance of X-ray spectral lines has been 
recently put into question by rigourous statistical argu- 
ments (|22). Therefore, we try to assess whether the line 
in Mrk 335 could be due to statistical fluctuations through 
Montecarlo simulations. The phenomenological model 
(power law H-broad em. line) used in the spectral analy- 
sis section, is taken as a "baseline" model and the Ax^ 
for adding an absorption line is measured to be 14.37 in 
the real spectrum. In this way, 10000 fake background- 
subtracted data sets have been obtained. Each of these 
spectra is fitted with the baseline model (power law + 
broad Gaussian line) and only then, a narrow absorption 
line with (t=50 eV is added to the fit, in order to measure 
the improvement in with respect to the value ob- 
tained by fitting the baseline model with no absorption. 
We obtain a Ax^ larger than in the baseline model in 263 
cases, yielding a significance for the absorption line of 
97.37 percent. 



3. SEARCH FOR ENERGY-SHIFTED NARROW 
LINES: SIMULATION PROCEDURE 



All narrow lines detected in the literature have been found 
in individual sources, as the one discussed for Mrk 335 
and many of them have been found marginal. Marginal 
detections could possibly arise due to random deviations 
in the spectra. Quantifying the significance of such devi- 
ations in a large number of X-ray spectra would provide 
an estimate of the robustness of the detections. To date, 
a systematic search for the presence of such features in a 
sample of objects has not been performed. An attempt to 
do that is presented in the following. A sample of archival 
PG quasars observed by XMM-Newton has been chosen. 
For completeness, the description of the X-ray proper- 
ties of this sample is reported by ( l23l) and (j^). In the 
present analysis instead, only the sources with more than 
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Figure 3. The plot shows the ratio of the computed 3-9 keV flux to a power-law with index T = 2.2 ( rest frame energy). 
Left: inflow with Tin — 20 Rg, rout — 10"^ Rg, — 0.2 Mq yr^^, neglecting radiation pressure. The absorption features 
are predominantly due to Fe XXVI Ka [around 6.5 keV] and K(3 [around 7.7 keV]. Right: inflow with Tin = 24 Rg, 
'''out = 48 Rg, $ = 0.3 Mq yr^^, accounting for radiation pressure due to electrons at 60 per cent of the Eddington 
limit. The absorption feature is due to FE XXVI Ka. The light histogram shows the observational data for Mrk 335, also 
normalised to a power-law flt 



800 counts above 5 keV are included in order to insure 
good statistics in the Fe K region. The list of the sources 
is shown in Tabled A blind search for positive and neg- 
ative flux deviations has been performed in the simulated 
spectra. This procedure is addressed to test for the pres- 
ence of a narrow line in each of the spectra in the sample. 
We shall distinguish between unshifted Fe K lines (i.e. 
emission lines in the range 6.4-7 keV) and those which 
we will call energy -shifted lines, i.e. any absorption line 
and emission lines shifted out of this range. 

• For each spectrum, 10000 spectra have been simu- 
lated with XSPEC assuming a baseline model with- 
out any narrow line, folding it through the same in- 
strumental response and adding Poisson noise. Such 
spectra have then been grouped according to the 
same criterion adopted for the real data set, i.e. 20 
counts per spectral bin. In this way, 10000 synthetic 
background-subtracted spectra were generated, with 
photon statistics corresponding to the exposure in 
the pn detector (see second column in TableQ. 

• Each of these spectra was fitted first with the base- 
line model. Then, each of them was fitted a second 
time with a model consisting of the baseline model 
plus a narrow line. The narrow line parameters were 
set as follows. Positive and negative deviations were 
allowed for the line flux. The width of the narrow 
line was fixed to 50 eV during the fitting and the en- 
ergy of the line was stepped in steps of 70 eV, cor- 
responding approximately to the instrumental reso- 
lution. To avoid any calibration uncertainties at the 
boundaries of the instrumental response, the line is 
searched across the energy range 2.5-9.5 keV. For 
each energy on this grid, the value of the Ax^ for 
adding the narrow line to the data was calculated 
with respect to the baseline model. When the mini- 
mum Ax^ was found, the values of the correspond- 
ing energy and line flux were recorded, in this way, 
the most significant narrow lines in the data were 



detected. 

• The presence of a narrow feature has been tested in 
the real spectrum applying the same grid of energy, 
so that the comparison is made consistently. 

• The procedure of fitting the spectrum with the base- 
line model and then adding a narrow line with the 
same grid of energies described above, has been re- 
peated for each simulated spectrum. In each of the 
fake data sets, the greatest improvement in for 
adding the narrow feature, Ax^(sim), is recorded, 
along with the energy and the flux of the line, pro- 
viding a list of 10000 detections. Then, the number 
of spectra where Ax^(sim) > Ax^(data) is counted. 
This quantifies the probability that an apparent fea- 
ture as significant as that in the real data could be the 
result of random noise. 



Two runs of simulations have been performed. The first 
run picks up the most significant line in each spectrum, 
with the adopted baseline model. The second run is per- 
formed for a limited number of spectra where an un- 
shifted emission Fe Ka line has been detected at >90% 
in the first run. In these cases, the Fe K line has been 
added to the baseline model and the procedure has been 
run again. Table |2] summarises the final results and the 
baseline model used for each spectrum. 

This table comprises the list of the 10 detections signifi- 
cant at >90% selected in the following way: 

• Energy-shifted features significant at >9Q% from 
the first run of simulations 

• Energy-shifted features significant at >90% from 
the second run of simulations (i.e. after including 
significant Fe Ka lines in the baseline model) 
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Table 1. List of the PG QSOs sample where the narrow 
lines blind search has been performed. The sources have 



Table 2. List of the shifted features detected at > 90% in 
each spectrum with the corresponding baseline model 



been fitted with a power law in the 2-10 keV range. 


Quasar 
- 


Exposure 

s 


Gamma 


Flux 
(10-1^ cgs) 


X^/d.o.f. 
- 


Mrk 335 


28401 


— - . -Xn iV> 

9 19^ 

-0.02 


11.81 


815/742 


IIIZW2 


10141 


, , Lfl fid 


6.78 


344/349 


IZw 1 


18189 




7.73 


529/450 


PG 0844+349 


9230 


-0.07 


4.69 


169/172 


PG 0947+396 


17392 


. _„_i_n fid 

1 90^ 


1.75 


197/217 


PG 0953+414 


10211 


— -1- n n 

"1 Ci'ii'-'-'-'^ 
^•"'-0.05 


2.92 


194/225 


PG 1048+342 


19801 


I 09 + 0. 06 


1.32 


196/197 


PG 1 1 14+445 


34902 




2.23 


363/386 


PG 1115+080 


37082 


'•"^-0.06 


0.26 


178/184 


PG 1202+281 


11448 


. „, + 0.05 
'•'^-0.05 


3.60 


216/241 


PG 1211+143 


46884 


1 7,+0.02 
-0.02 


3.06 


741/638 


PG 1407+265 


46124 


24C+O.O2 
•'■^■'-0.02 


1.34 


465/490 


PG 1415+451 


20569 


1 QO+0.07 
'■"-0.08 


1.09 


151/159 


PG 1425+267 


29394 


, 40+0.04 
l-™-0.04 


1.65 


337/312 


PG 1427+480 


30763 


^•'°-0.06 


1.04 


246/251 


Mrk 478 


18037 


J 19 + 0.06 
-0.06 


1.84 


189/202 


PG 1448+273 


17840 


9^+0.06 
■"-0.06 


2.09 


239/245 


Mrk 841(*) 


7609 


, q,+0.03 

-0.03 


14.80 


468/492 


PG 1512+370 


13503 


, si+0.06 
'■°'-0.06 


1.86 


226/216 


PG 1634+716 


12207 


2 ,,+0.06 
■^"^'' — 0.06 


0.74 


267/207 


UGC 11763 
Mrk 304 


23579 
6945 


1 63+°°* 

f, 00+0. 06 


3.68 
3.41 


361/423 
224/115 



(*)3 spectra in the archive 



In total, 10 detections out of 24 spectra have been found 
at a significance higher than 90%. 

To estimate the probability to detect 10 features by ran- 
dom chance, it is assumed a null hypothesis in which the 
sample comprises 24 featureless spectra. The probability 
that this hypothesis can be rejected in the present data is 
then calculated assuming the binomial distribution as a 
probability distribution: 



Quasar 


Baseline 


Energy (*) 


Intensity (*) 


Significance 






(keV) 


(Phs^i cm"^) 




Mrk 335 


plaw+broad gau 


5.93 


-7.422e-06 


95.82% 


111ZW2 


pi aw 


9.01 


-5.870e-06 


90.93% 


IZw 1 


plaw 


3.20 


-1.187e-05 


96.34% 


PG 1211+143 


plaw+zwabs 


7.61 


-2.9541e-06 


100% 


PG 1211+143 


plaw+zwabs+2 zgau 


2.92 


-7.411e-06 


99.96% 


PG 1407+265 


plaw 


7.89 


3.794e-06 


92.2% 


Mrk 841 (I) 


plaw 


6.28 


1.504e-05 


98.26% 


PG 1634+716 


plaw 


4.11 


-1.109e-5 


96.99% 


UGC 11763 


plaw 


6.28 


3.781e-06 


94.21% 


Mrk 304 


plaw+zwabs +zgau 


4.46 


9.986e-06 


100% 



is rejected. From statistical considerations, it is reason- 
able saying that out of 10 detections some are real. From 
a more speculative point of view, let us considering the 
number of successful events characterised by the aver- 
age cumulative probability Pc=0.5. At this point, it is 
as likely to have more than x false detections, as to have 
less than x false detections, because obviously they have 
the same cumulative probabihty. For 24 number of tri- 
als, the probabihty calculation shows that such number 
is between 2 and 1 . So, one could draw an approximate 
conclusion by saying that the majority of the 10 detec- 
tions are real, at a level of confidence of 90% and 2-3 of 
them are false. 



4. SUMMARY OF RESULTS 



n! 



x\{n — a;)! 



(2) 



In the present case, the parameters of the distribution cor- 
respond to n=24 (number of spectra), x=10 (number of 
successes i.e. detections at >90%), p=0.l (i.e. the prob- 
ability to have a detection at more than 90% in the case 
of the null hypothesis). The cumulative probability for a 
given number of random detections defines the probabil- 
ity to have up to that number of random detections in the 
sample and it is defined as: 

X 

fc=0 

The probability to find 10 or more random detections in 
the sample is calculated to be 1-Pc=5.25 x 10^^. This 
number is very small, implying that it is very unhkely 
that 10 detections in a sample of 24 spectra are all due 
to random noise. Therefore, the possibility that none of 
them is real can be ruled out, i.e. the null hypothesis 



A narrow absorption line has been detected at ~5.9 keV 
with a significance of ~97% in the EPIC pn spectrum of 
Mrk 335; if interpreted as Fe XXVI Ka and if the ef- 
fect of the gravitational field is neglected, the observed 
redshift of the line corresponds to a receding velocity of 
50000 km s~^ in the absorbing gas. The comparison to a 
physical inflow model shows that the line is consistent 
with being produced in a discontinuous flow of mate- 
rial dragged in high velocity motion towards the nucleus, 
rather than in a spherical flow. Arguably, the fact that 
the line is not smeared nor broad is a strong indication 
against the hypothesis of matter in orbit at a few gravita- 
tional radii, as suggested for other similar cases. 

The blind search for energy shifted features carried out in 
the sample of PG quasars provided an encouraging result 
on the statistical significance of the narrow lines, in gen- 
eral. The majority of the detections are in fact believed 
to be real in the XMM-Newton data. This prehnainary re- 
sult should be investigated using a much larger sample of 
spectra. 



6 



ACKNOWLEDGMENTS 



The AstroGroup at Imperial College London is acknowl- 
edged for financial support. A.L.L. is grateful to Gio- 
vanni Miniutti for many stimulating discussions during 
this conference. 



REFERENCES 

[I] Nandra, K., George, 1. M., Mushotzky, R. E, Turner, 
T. J., and Yaqoob, T. ASCA Observations of Seyfert 1 
Galaxies. II. Relativistic Iron K alpha Emission. ApJ, 
477:602-1-, March 1997. 

[2] Nandra, K., George, 1. M., Mushotzky, R. R, Turner, 
T. J., and Yaqoob, T. The Properties of the Relativis- 
tic Iron K-Line in NGC 3516. ApJ, 523:L17-L20, 
September 1999. 

[3] Yaqoob, T. and Serlemitsos, P. Iron K Features in 
the Quasar E1821H-643: Evidence for Gravitationally 
Redshifted Absorption? ApJ, 623:112-122, April 
2005. 

[4] Longinotti, A. L., Cappi, M., Nandra, K., Dadina, 
M., and Pellegrini, S. The complex FeK Une of 
the Narrow-Line Seyfert 1 galaxy IRAS 13349H-2438. 
A&A, 410:471-479, November 2003. 

[5] Ruszkowski, M. and Fabian, A. C. On the influence 
of resonant absorption on the iron emission-line pro- 
files from accreting black holes. MNRAS, 315:223- 
228, June 2000. 

[6] Matt, G., Porquet, D., Bianchi, S., Falocco, S., 
Maiolino, R., Reeves, J. N., and Zappacosta, L. A 
changing inner radius in the accretion disc of Q0056- 
363? A&A, 435:857-861, June 2005. 

[7] Dadina, M., Cappi, M., Malaguti, G., Ponti, G., and 
De Rosa, A. X-ray absorption lines suggest matter in- 
falling onto the central black-hole of mrk 509. astro- 

ph/0506697, 2005. 

[8] Reeves, J. N. et al. Evidence for gravitational infall 
of matter onto the super- massive black hole in the 
quasar pg 1211+143? 2005. 

[9] Turner, T. J., Mushotzky, R. F, Yaqoob, T, George, 
1. M., Snowden, S. L., Netzer, H., Kraemer, S. B., 
Nandra, K., and Chelouche, D. Narrow Components 
witiiin the Fe Ka Profile of NGC 3516: Evidence of 
the Importance of General Relativistic Effects? ApJ, 
574:L123-L127, August 2002. 

[10] Turner, T. J., Kraemer, S. B., and Reeves, J. N. 
Transient Relativistically Shifted Lines as a Probe of 
Black Hole Systems. ApJ, 603:62-66, March 2004. 

[II] Guainazzi, M. The history of the iron K Une profile 
in the Piccinotti AGN ESQ 198-G24. A&A, 401:903- 
910, April 2003. 



[12] Yaqoob, T, George, 1. M., Kallman, T. R., Padman- 
abhan, U., Weaver, K. A., and Turner, T. J. Fe XXV 
and Fe XXVI Diagnostics of the Black Hole and 
Accretion Disk in Active Galaxies: Chandra Time- 
resolved Grating Spectroscopy of NGC 7314. ApJ, 
596:85-104, October 2003. 

[13] Porquet, D., Reeves, J. N., Uttley, P., and Turner, 
T. J. XMM-Newton observation of the Seyfert 1.8 
ESQ 1 13-GOlO: Discovery of a highly redshifted iron 
line at 5.4 keV. A&A, All: 101-105, November 2004. 

[14] Gallo, L. C, Fabian, A. C, BoUer, T, and Pietsch, 
W. A possible line-like emission feature at 8 keV in 
the Seyfert 1.2 UGC 3973. MNRAS, pages 732-+, 
August 2005. 

[15] Nayakshin, S. and Kazanas, D. Narrow Moving Fe 
Ka Lines from Magnetic Flares in Active Galactic 
Nuclei. ApJ, 553:L141-L144, June 2001. 

[16] Dovciak, M., Bianchi, S., Guainazzi, M., Karas, 
v., and Matt, G. Relativistic spectral features 
from X-ray-iUuminated spots and the measure of the 
black hole mass in active galactic nuclei. MNRAS, 
350:745-755, May 2004. 

[17] Gondoin, P., Orr, A., Lumb, D., and Santos-Lleo, 
M. XMM-Newton observations of the Seyfert 1 
galaxy Mrk 335. A&A, 388:74-87, June 2002. 

[18] Fabian, A. C, Rees, M. J., Stella, L., and White, 

N. E. X-ray fluorescence from the inner disc in 
Cygnus X-1. MNRAS, 238:729-736, May 1989. 

[19] Nayakshin, S. and Kallman, T. R. Accretion Disk 
Models and Their X-Ray Reflection Signatures. I. Lo- 
cal Spectra. ApJ, 546:406-418, January 2001. 

[20] Sim, S. A. ModelUng the X-ray spectra of high- 
velocity outflows from quasars. MNRAS, 356:531- 
544, January 2005. 

[21] Gierlihski, M. and Done, C. Is the soft excess in 
active galactic nuclei real? MNRAS, 349:L7-L11, 
March 2004. 

[22] Protassov, R., van Dyk, D. A., Connors, A., 
Kashyap, V. L., and Siemiginowska, A. Statis- 
tics, Handle with Care: Detecting Multiple Model 
Components with the Likelihood Ratio Test. ApJ, 
571:545-559, May 2002. 

[23] Piconcelli, E., Jimenez-Bailon, E., Guainazzi, M., 
Schartel, N., Rodriguez-Pascual, P. M., and Santos- 
Lleo, M. The XMM-Newton view of PG quasars. I. 
X-ray continuum and absorption. A&A, 432:15-30, 
March 2005. 

[24] Jimenez-Bailon, E., Piconcelli, E., Guainazzi, M., 
Schartel, N., Rodriguez-Pascual, P. M., and Santos- 
Lleo, M. The XMM-Newton view of PG quasars. 
II. Properties of the Fe Ka line. A&A, 435:449-457, 
May 2005. 



